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part of the matter-induced gravitational couplings of the neutrinos under such conditions, which 
arise at order g 2 /M^ Vl and determine their contribution to the neutrino dispersion relation in 
the presence of a gravitational potential cxt . These new contributions vanish for the muon and 
tau neutrinos. For electron neutrinos with momentum K, they are of the order of the usual 
Wolfenstein term times the factor (K 2 /M|^)</> ext , for high energy neutrinos. In environments 
where the gravitational potential is substantial, such as those in the vicinity of Active Galactic 
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be taken into account in the analysis of possible violations of the Equivalence Principle in the 
neutrino sector, in experimental settings involving high energy neutrinos traveling through a 
matter background. 



Jose F. Nieves 

Laboratory of Theoretical Physics 
Department of Physics, P.O. Box 23343 
University of Puerto Rico, Rio Piedras 



Puerto Rico 00931-3343 



Palash B. Pal 



Theory Group, Saha Institute of Nuclear Physics, 
1 /AF Bidhan-Nagar, Calcutta 700064, India 



Abstract 



1 Introduction 



In the vacuum, a photon that propagates with momentum K in a constant gravitational potential 
cj) ext has a dispersion relation given by 




(1.1) 
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In fact, this formula holds not just for photons, but for any massless particle. In particular, it also 
holds for neutrinos in the limit in which their mass can be neglected. This result is a consequence of 
the universality of the gravitational interactions that is embodied in the Principle of Equivalence, 
according to which gravity couples to matter particles in a universal way; i.e., the graviton field 
couples to the stress-energy tensor of the particles with a coupling constant that is the same for all 
particle species. 

Since the gravitational contribution to the neutrino dispersion relation is the same for all the 
neutrino flavors, it is not relevant in the context of neutrino oscillation experiments. That contri- 
bution yields a common factor in the evolution equation and therefore it drops out of the formulas 
for the flavor transition amplitudes. 

On the other hand, it has been pointed out by Gasperini and by Halprin and Leung §, that 
a small violation of the equivalence principle, manifest as a difference in the couplings of the various 
neutrinos to gravity, would have observable consequences in various neutrino experiments. By the 
same token, precise measurements of various observable quantities in these experiments can place 
significant constraints on any deviations of the Equivalence Principle in the neutrino sector^]. 

In our earlier work [Q, || we considered neutrinos traveling through a background of normal 
matter, made of electrons and nucleons. By calculating the weak interaction corrections to the 
stress energy tensor of the neutrinos in the background, we showed that the lowest order gravita- 
tional couplings of the neutrinos are modified in a non-universal way. In summary, the magnitude 
of the background-induced contribution to the neutrino gravitational couplings was calculated to 
order g 2 /M^r, the corrections to the neutrino dispersion relations in the presence of a gravita- 
tional potential were determined, and some of their possible phenomenological consequences were 
considered Q. 

The corrections to the neutrino effective potential determined in Ref. 0] are independent of the 
neutrino momentum because, to order g 2 /Myy, the weak interactions are local. Our purpose in the 
present work is to extend the calculation of the stress-energy tensor of the neutrinos, to include the 
order g 2 /M^ terms. These corrections arise by retaining the momentum-dependent terms of the 
VK-propagator in the relevant one- loop diagrams for the neutrino gravitational vertex in matter. 
The procedure is analogous to that employed in Ref. [0, ^] to determine the non-local corrections 
to the neutrino effective potentials in matter in the absence of any external field. 

The distinguishing feature of the 0(g 2 /M&r) corrections that we determine in the present work 
is that they induce momentum-dependent terms in the effective potential of the neutrinos in the 
presence of a gravitational potential, which mimic the features of a non-universal gravitational 
coupling of the neutrinos at a fundamental level. We point out that our calculations, and the 
results based on them, do not depend on any physical assumption beyond those required by the 
standard model of particle interactions and the linearized theory of gravity, including the question 
of whether or not the neutrinos have a non-zero mass. Hence, the effects that we will consider are 
present at some level and it is conceivable that they are detectable in some favorable situations 
involving strong gravitational fields, such as those that exist in the vicinity of active galactic nuclei. 
In any case, these corrections in principle must be taken into account in the analysis of the tests 
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and constraints of possible violations of the Equivalence Principle in the neutrino sector, related to 
experimental observations involving neutrinos traveling through a background, as opposed to the 
vacuum. 

The paper is organized as follows. In Sec. ^, we present the diagrams to be calculated and the 
Feynman rules needed for their calculation. We write down the expressions for vertex diagrams 



with internal Z-boson and VF-boson lines in Sec. |3.l| and Sec. |3.2| respectively. This sets up the 
stage for Sec. [|, which contains most of the results of this paper. In this section, we start by 
setting up the general formalism which relates the neutrino gravitational vertex function to the 
neutrino self-energy in an external gravitational field. Then, in Sec. 4.2 , we show that although the 
general vertex can depend on the unphysical parameter £ appearing in the W and Z propagators, 
the dispersion relation is independent of it. The diagrams with internal Z-bosons lines do not 
contribute to the dispersion relations at 0(g 2 jMy/). The 0(g 2 /Myy) contributions of the different 
W-mediated diagrams to the self energy are calculated in Sec. and they are collected together 
in Sec. || to find the neutrino dispersion relation to 0(g 2 /My^). We conclude with some comments 
in Section H. 



2 Feynman rules 

The diagrams that we must calculate in this paper are shown in Figs. [I] and |2[ In principle, we 
must consider also the diagrams in which any number of internal vector bosons is replaced by 
the corresponding unphysical Higgs bosons. In general, only by considering them we can attempt 
to prove that their contribution cancels the one that arises from the gauge-dependent piece of 
the gauge boson propagator. However, for simplicity, we will carry out the entire calculation in 
the limit m e — ► 0, in which case the diagrams involving the unphysical Higgs particles do not 
contribute. Hence, we do not consider them any further. Our task is therefore twofold: to calculate 
the 0(1/ Mg) contribution terms of the diagrams shown in Figs. [I] and ^ where B is either the W 
or the Z-boson, and to show that the contribution from the gauge-dependent terms of the vector 
boson propagators are zero in the limit m e — > 0. 

The propagators and vertices that we need to perform the computation are as follows. The 
thermal propagator for a fermion of mass mt is given by 



iS f(p) = {j> + m f) 

where 



p 2 — m f + ie 



1 ——-2tt5{p -m f )rj f (p) 



(2.1) 



Vf( P )= + ff-P'? (2.2) 

with (3 being the inverse temperature, [if the chemical potential and the velocity four- vector of 
the medium. 
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Vi(k') Ui(k) Ui(k') ui(k) u^k') Vi(k) u^k') Vi{k) 

(A) (B) (C) (D) 

Figure 1: Z-exchange diagrams for the one-loop contribution to the gravitational vertex of any neutrino 
flavor (i = e, /x,r) in a background of electrons and nucleons. The braided line represents the graviton. 




Figure 2: VF-exchange diagrams for the one-loop contribution to the v e gravitational vertex in a 
background of electrons. 
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The propagators for a gauge boson B = W,Z are given in the £ gauge as 

The contributions of order to the neutrino effective gravitational vertex were calculated in 

Refs. ||[§. Here we are interested in the terms next to the leading ones, of order 1/Mg. For that 
purpose, it is sufficient to approximate the boson propagators that appear in the loop integrals by 

Aj^fo) = Ig + Ag^ + Ag^), (2.4) 



where 



4g» = (2-5) 

The remaining terms would yield contributions of order 1/Mg to the effective vertex, which we are 
not considering. 

B v {k') 




™M%a' Xpilv - in (C xpiiu (k, k>) + C^(k, k')) 



Figure 3: Feynman rule for the gravitational vertex coupling of a gauge boson B = W,Z with the 
graviton. 

We also need the Feynman rules for the graviton couplings of the vector bosons, including 
the dominant sub-leading terms that were not considered in our earlier work, as well as the gauge- 
dependent terms. To obtain them we need to start from the free Lagrangian, which for the VF-bosons 
is 

=A/^{-^r<^Wp^ , (2.6) 

where g = det(gw). We now write 

9iw = Vtiv + 2nh pi/ , (2.7) 
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where k is related to the Newton's constant by 

K = V8wG, (2.8) 

and hav is the graviton field. Keeping up to first order terms in k, we obtain the VFPF-graviton 
coupling in the form 

^ ww ) = - K h^(x) (f Xp (x) + Tj§(xj) , (2.9) 

where 

Tx P {x) = \vx P W; v W^ + 2W^W%-M^a! XpiiV W* v W» 

TxJ(x) = s [nx P (d ■ w*)(d • w) - (d x w; + d p w* x )d -w-d- w*(d x w p + d p w x )} ,(2.10) 

using the shorthand 

a 'x P pu = VXpVfj,u ~ (VXuVpu + VXvVph) ■ (2-11) 

The Feynman rule for the VFW-graviton coupling can then be written as in Fig. |3], with a Xppu as 
in Eq. (pUP , and 

Cxpfiv{k, k ) = ^xpiVp.fk • k — k^k^) — Tj pu {kxk p + k x k p ) 
+ K{ri Xp k' p + rj pp k' x ) + k'^(ri Xu k p + r] pv k x ) 
- k ■ k'(rj Xfl r] pv + r] Xl/ T] pfi ) 

^Ap/^(^> = £ {vxpkpk'^ — k' x k p r] pu — k l p k p r\ Xv — k x k' v rj pp- — k p k' v rj Xp ^ . (2.12) 

The rule for the couplings of the Z boson are the same. 

The other couplings necessary for the present calculation have been deduced in Ref. [||. For the 
gravitational coupling of a fermion of mass nif, the Feynman rule involves the factor —inV Xp (p,p'), 
with 

V\ p (p,p' ) = j lx{p + p') P + l P {p + p')x -^ + p"-2m f , (2.13) 

where p and p' denote the momentum of the incoming and the outgoing fermion, respectively. 
There is also a fermion vertex involving both a graviton and a gauge boson line, and its Feynman 
rule is given in Fig. with the definition 

d.Xpa/3 = VXpVaP - ^{VXaVp/3 + VpaVXp) ■ (2-14) 
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Figure 4: Feynman rule for the graviton coupling involving the Z-boson and a fermion. a\ pa p 
has been defined in Eq. ( 2.14| ) and, in the standard model, X u = —Y v = 1/2 for a neutrino, 



and -X e = X p = \ - 2sm 2 6 w , X n = -\, Y e = \, Y n = -Y p = \g A , with g A = 1.26 being 
the renormalization constant of the axial-vector current of the nucleon. The analogous couplings 
involving the W boson are obtained by replacing 2 cAe w ^ V"2' ^ 3 anc ^ ^ ^ — 3" 

3 Calculation of the vertex function diagrams 
3.1 The Z-mediated diagrams 

The diagrams, which are shown in Fig. [I], are the same as in Ref. [Q]. We use the same notation and 
conventions given in that reference for the weak-neutral-current couplings in the Lagrangian [e.g., 
Eqs. (2.31)-(2.35) there], summarized in Fig. || here. The Z-propagator is external to the loop, so 
all the algebraic manipulations of the loop integral are performed exactly as in Ref. pj . As in that 
reference, we define 

b^=J2x f b f , (3.1) 
/ 

where the sum is over the fermions species in the background, andQ 

b f = AV2G F J ^5{p 2 - m 2 } )r, f { P )p-v 

= \/2G F (n f - nj) , (3.2) 

with rijj being the number densities of the particles and antiparticles respectively. We denote 
by k, k' the momentum of the initial and final neutrino respectively, and the momentum of the 
outgoing graviton is 

q = k-k' . (3.3) 

For Fig. |l]A we then obtain 

rr = (7« + «V-(i-«^ ( A g> _ b ^ fVa) L , (3 . 4) 

1 We take the opportunity to correct a typographical error in Eq. (3.20) of Ref. H 



7 



where 



IV Xpa 



9z 

Ml 2 ? 



d 3 p 



2E f {2Trf 



q 2 -2p r q 



+ {q^ -q) 



-Yf(ff + h) 



N ( xJ a (Pf,q) 
q 2 — 2pf • q 



(q -> -q) 



(3.5) 



with ff and fj denoting the momentum distribution functions for particles and antiparticles and 



N \pa(p>i) = ( 2 P-q)\[ 2 PpPa- (p a q P + q a p P ) + (p ■ q)r/ a p} + (\ <^ p) 



N { xpa(p,<l) = (2p - q)\ie pa ^V a + (A <-> p) . 



(3.6) 



In Fig. [TjB, the Z line attached to the fermion loop carries zero momentum, in which case the 
C\pfj, u and terms of the ZZ graviton couplings do not contribute. Hence the contribution 

from this diagram is given by 



Xp 



Xpa/3 7 + 



q 2 r-(l-Qq c 
M 2 



(3.7) 



where a' x a has been defined in Eq. (|2.11 



In Fig. [I]C, the Z-boson line has zero momentum flowing through it. As a consequence, the 
0(M^ 4 ) terms of the Z-propagator do not contribute and the amplitude for this diagram is the 
same as obtained in Refs. [|], ||, viz., 



n(lC) 



b (Z) 



aXpa/31 



where a\ pa p was defined in Eq. ( p.!4| ). Finally, for Fig. [l[D, we get 



Xp 



b (Z) 



(** + - c 1 - ^° 

OApo/3 7 H 



Mf 



(3i 



(3.9) 



In this way, summing up all the diagrams of Fig. [T], we obtain the total contribution from the 
Z-diagrams to the effective vertex 



v(l) 

Xp 



+ 



fr - (1 - e)g c 

Mf 



VXpVaf3 + 7;(VXaVpf3 + VpaVxp) 



(3.10) 



where the first term on the right hand side represents the leading order contribution, of order 
(1/Ml), whose detailed evaluation was the subject of Ref. [|]. It was shown in that reference that, 
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for a static and constant gravitational potential, which involves taking the q — > of the vertex 

(Z) 

function in the proper way, the contribution from the A^a term to the neutrino self-energy is a 
well defined quantity. Borrowing that result here, it follows that the remaining terms in Eq. ( |3.10| ), 

(Z) 

which involve the same quantity , yield a contribution to the neutrino self-energy that depends 
on the derivatives of the gravitational field. Therefore, in the presence of a static and homogeneous 
potential, the 0(1/ M^) contributions to the neutrino self-energy vanish. Since these are the only 
diagrams that contribute to the tv T dispersion relation if there is no fi or r in the background, as 
we assume, only the dispersion relation of the v e is modified by the 0(1/Mjj) terms. These arise 
from the W diagrams that we consider next. 

3.2 The IV-mediated diagrams 

We now discuss the diagrams involving the VF-bosons, shown in Fig. [2l which are relevant only for 
the electron neutrino. In what follows, unless we explicitly mention it otherwise, we are referring 
only to this type of neutrino. 

3.2.1 Diagram §A 

Let us start with Fig. §A. It gives 

= -\j j^l a LiS e (p f )V Xp (p, P f )iS e (p)^LA WaP (p - k) , (3.11) 

where 

p'=p-q. (3.12) 
By substituting the W propagator in Eq. ( |3.11 ) three terms are produced. The 0(1/M^) term 



has already been calculated in Ref. Qj. Separating it out, we can write the dispersive part of T'^ 4 ^ 
in the form 

rr = Air-M^+^+<), (3.i3) 

where , which contains all the 0(1/ My^) contributions that are £- independent, is given by 



(3.14) 



while H[f is given by the same expression as Eq. (|3.14| ) but with the replacement 



A$*-A<00. (3.15) 
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3.2.2 Diagram gB 

Proceeding in similar fashion for this diagram, we borrow the leading order result derived in Ref. 
and write 



r' 



(2B) 
Ap 



1 



L + <>+<> 



(3.16) 



( B) ( B) 

where G\ p and H^ p contain the non-leading terms, separated according to whether they are £- 
independent or not, respectively. These non-leading terms can come either from the non-leading 
contribution in the PFiy-graviton coupling or in the TV-propagator. Accordingly, we get the fol- 
lowing terms 



G 



Xp 



a 



d 4 p 
(2^)3 

-^C XpaP {p -k', P -k)+ ^ai w A^(p — k) + TTa'^A^lp - k>) 



H 



(B) 
Xp 



1 w 

d 4 p 



(2ny 



5(p 2 ) Ve (ph a L^L 



1 



Cf^p -k',p-k)+ fo'^lp - k) + tT^A^Up - k') 



(0 



(3.17) 



3.2.3 Diagrams |C and D 

Similarly, we define the non-leading contributions by writing 



r (2C) 
1 Xp 



^ + j(7a^7p + 7 p ^7a) 



(3.18) 
(3.19) 



For these diagrams, the non-leading terms come only from the TV propagator. Thus we obtain 

-,2 



G 



G 



(c) 

Ap 

(D) 
Xp 



9 

y fl Ap/3p 



'(^Xpap 



^6(p 2 )rj e (p) la L^LA^ a (p - k') 
d 4 p 



(3.20) 
(3.21) 



The corresponding formulas for H^' D ' are given by the same expressions, with the replacement 
indicated in Eq. ( 3,15| ). 
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4 Self-energy in a constant potential 



4.1 General formalism 

In the presence of a static and homogeneous gravitational potential, the neutrino self-energy, in- 
cluding the 0(g 2 /Myy) corrections, is given by 

S = S mat + 4 2) + 4 4) + sg ) (4.1) 

where S mat stands for the purely weak self-energy terms which do not depend on the gravitational 
potential, including the 0{g 2 /M^) corrections, while Ytq is the quantity calculated in Ref. 
[defined in Eq. (4.8) there]. In analogy with the expression given in Eq. (4.9) of that reference, we 
can write here 

4V,i?) = cxt (2^V-^)G Ap (O,Q^O), 

= cxt (2v x vP-ri x P)H Xp (O,Q^O), (4.2) 

where G\ p (Q, Q) and H\ p (£t, Q) denote the sum of the various gauge-independent and gauge- 

(X) (X) 

dependent terms G\ p and H\ ' respectively, whose integral expressions have been given in Eqs. 



( fTlll ), ( fmD , ( |3720| ) and ( |3T2l] ) . As we already argued at the end of Sec. |0|, the 0(g 2 /Mf) terms 



from the Z-mediated diagrams do not contribute to the self-energy of the neutrinos, and we have 



to calculate the contributions only from the VF-mediated diagrams in Eq. (4.2). In Eq. (4.2) we 

have indicated the dependence of the various quantities^ on the kinematic variables of the problem. 
These variables are defined by writing 

= (uj,K) 

<f = (H,Q) (4.3) 

in the rest frame of the medium where 

^ = (1,0). (4.4) 

In what follows, we indicate the dependence of the various functions on these variables when it is 
necessary, but we omit it otherwise. 

The chirality of the neutrino interactions dictate that £ has the decomposition 

S = (a^ + 6^)L, (4.5) 

where a and b are in general functions of co and K. It is convenient to split a and b as 

a - n W 2) +n (4) + 
a — Omat ~r Qq + a G + O g 

b = 6 ma t + 4 2) + b G + bi G> ( 4 - 6 ) 



2 Although we do not show it explicitly, G' x ' depends on uj and K as well. 
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corresponding to the decomposition given in Eq. (|4.l|). Eq. (|4,5| ) implies that a G and b G can be 
determined by means of the formulas 



,( 4 ) 



4 4) = - „ 2 "" , (4-7) 



if 2 

^(4) _ fc 2^(4) 
K 2 



and similarly 



a G 



where 



B^(uj,K) = ext (2i;V-^)^G ! Ap (O ) Q^O) 



and similarly 



A^(lo,K) = ^\2v x vP-r 1 x P)\^Gx p (^Q^Q), (4.£ 



B^,*) = cxt (2 V V- r^Tr^^O, Q ^ 0) , 

= (f^V-r/^Tr^^Q^O). (4.10) 



As discussed in Ref. Q, the dispersion relations are obtained by solving 



uok = K + 



1 — a(tJK, K) 



OK = K - K 7 K - A 1, (4.11) 
1 - a(-uj K ,K) 

for the neutrinos and antineutrinos, respectively. In the context of our perturbative approach, the 
solutions are given approximately by 

u K = K + [1 + a(K, K)]b(K, K) 

u K = K —[1 + a{—K, K)]b(—K, K) . (4.12) 

Using the result for the dispersion relation obtained in Ref. Q [given in Eq. (4.32) of that 
reference], it is easy to see from Eq. (4/7) that in the present case 
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Here b(? is the contribution calculated in Ref. @ [denoted by simply be and given by the formula in 
Eq. (4.30) there], while 

&mat and 6 ma t are the usual contributions independent of the gravitational 
potential, for neutrinos and antineutrinos respectively, including the 0(g 2 /Myy) terms. These can 
be obtained from the references cited [0, H|, and we will quote explicit formulas for our cases of 
interest in Sec. pi Thus, what remains is simply to calculate the coefficient by means of Eq. 
( |4.9| ), using the integral formulas for the various functions G\ obtained in the previous section. 
However, before we proceed to do that, we will now show that the quantity B^(±K, K), which is 
calculated in similar fashion, is in fact zero so that we finally obtain the gauge invariant result 



K + 
K + 



iext 



K + Q. + 0& mat + # + 



b ^ K + (i + <f xt )& mat - bg> + 



BW(K, K) 
K 

B^(-K,K) 
K 



(4.14) 



4.2 Gauge independence of dispersion relations 

We calculate B^ using Eq. ( [4.10|) . We denote by B^ the contribution from each term that 
arises in that formula. Carrying out the trace operation first, a straightforward calculation yields 
the expressions 



Bi A \ 



B<f\u,K) 



Bf\u,K) + B^>(u J ,K) 



(D), 



icxt 



icxt 



i cxt 



g 2 i 
24, 



k 2 (2v x vP - rfo) [X« + Y Xp (0, Q 0) 



k 2 (2v x v p -7] xp ) 



a Xpa(3 



X W"P + X 



(3) 
Xp 



k 2 (2v x v p - r] X P)a Xpal3 X^ 



where 



X 
X 

t 

X 



(i) 
Xp 

(2) 
a/3 

(3) 



d 4 p 



(2vr) 3 
d 4 p 



S(p 2 )r] e (p) 



- } r]Xp{k • p) + 2p x pp 



(2vr) 3 
d 4 p 



HP )Ve(p) [2paPf3 -Pakp- Ppk 



/dp 
■7^-r^S(p 2 )r] e (p) [rjxpP ■ k + 4p x p p - 2p x k p - 2p p k x ] 



(2vr) 3 
d 4 p 



Hp 2 )Ve(p) 



JXp 



2p ■ q 



with 



Y Xp (n, Q) 

ix P = j(2p- q) x [(2k-p-q ■ k)p p - (k ■ p)q p + (p • q)k p ] + (A 



P) 



(4.15) 



(4.16) 



(4.17) 
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(1 2) 

The terms ' can be expressed very simply in terms of integrals over the electron distribution 
functions while Y\ p (0, Q — > 0) can be evaluated following the techniques used to evaluate similar 
integrals in Ref. Q. However, for our purposes here, the explicit results for them are not needed. 
For consistency, we just note that, despite appearances, (2uV - 7] X P)Y Xp (0, Q -> 0) is well-defined 
(see Appendix [A]). Since each expression in Eq. ( 4.15[ ) contains an explicit factor of fc 2 ,we then 
obtain 

Bt(±K,K) = 0, (4.18) 



which completes the proof of Eq. ( 4.14 ). 
4.3 Calculation of 

As a byproduct of the gauge invariance proof given above, it follows that the longitudinal part of 
does not contribute to the dispersion relations, which involve B^(±K, K). Therefore, to 
obtain the dispersion relations, we can use 



k 2 = u? - K 2 = (4.19) 



and substitute 



in the formulas for the G^K Corresponding to each there is a term in and according 

to Eq. ( |4.9|) we have 

B^(K,K) = oxt ( C (x) (0, Q -> 0) ) , (4.21) 

y X=A,B,C,D 



where 



G W (O, Q) = (2v x v? - rj^Tr^Gif . (4.22) 



We consider the computation of each of the G^ x \i}, Q) separately. 
4.3.1 Computation of G^ 



In the term containing r) e (p') in the integral given in Eq. ( |3.14 ), we make the change of variables p 



p + q. Then using the replacement indicated in Eq. ( 4.20[ ) for the W propagator, and remembering 
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that we are setting m e 



where 



0, we obtain 
g 2 f d 4 p 



2M^ 



(2vr)- 



(p-k) 2 
I 2 — 2p ■ q 



l\2 



+ 



q 2 + 2p ■ q 



^\ P {q) 



-Tr 



-Tr 



fflV\ p (p,p- q)(j/>- I 
fflV\ p (p,p- q)(j/>- fH)j5 



(4.23) 

(4.24) 
(4.25) 



In writing these forms, we have used the cyclicity of trace and the contraction identity 7^7" = —2$. 
We have also used the identity Cj^C -1 = — 7J to invert the order of the gamma matrices inside 
the trace in the second term in Eq. ( 4.23f) , as well as the property V\ p (p,p') = V\ p (p',p). The trace 
operation can be carried out straightforwardly. The ^t\ p terms does not contribute to the final 
result. As for the *dK\ p terms, we carry out the integral over p° with the help of the 5-function and 
substitute the resulting expression in Eq. ( gjgp . This gives 



d 3 p 



+ 



k 



r> 



(2tt) 3 2£ 
2k' -p 



2k -p 



,q)fe + 



k 2 + 2k ■ p 



q 2 + 2p ■ q 



Jt(p,-q)fe + 



2p ■ q 

k' 2 + 2k' ■ p 



--k-p 



k 2 + k' 2 



2p ■ q 
1 



q 2 + 2p ■ q 

■^(-P,-Q)fe 



-p,q)fs 



2p ■ (k + k')\ f e + -k-p 



k 2 + k' 2 + 2p-(k + k')] f, 



(4.26) 



where 



JZ{p, q) 



+ 



4(p • v) 2 + 4p ■ v q ■ v — (q ■ v) 2 k ■ p 

— 2p-v + q- vk-vp-q+ 2p ■ v — q ■ v p ■ v k ■ q 



(4.27) 



To evaluate G^ A \0,Q — > 0) we proceed as follows. We carry out the integration in the rest 
frame of the background. In that frame Eqs. (|4.4j ) and (|4.3j ) hold, and similarly we write 



(S,V) 



(4.28) 



where £ = |"P| because of the (5-function that appears in Eq. ( 4.23j ). In addition we use the fact 
that k^ satisfies Eq. ( 4.1SQ . We next evaluate Eq. ( 4.26| ) in the static limit O = and take Q — > 
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afterwards. The terms without any denominator can be easily evaluated in this limit by putting 
k = k' . The results can be conveniently expressed by defining the integrals 



7W"Pa* fe(£) + (-l) r h(£) 



(2tt) 3 2£ 

so that the contribution of these terms to G^ A \0, Q —* 0) can be written as 



(4.29) 



Ik^k I^ii . 



(4.30) 



It is now useful to introduce the scalar integrals, for any fermion /, 

d 3 V 



/(/) 



1 

4^2 



(2tt) 3 2£ £ 
d£ £ n+1 



f f (f) + (-l) n //(£) 

/)(£) + (-!)"/#) 



(4.31) 



for which we give explicit formulas in Sec. 0. Further simplification is then obtained by noting that 
the integrals in Eq. (4.29) can be determined in terms of the J„ ; e.g., 



J (e) W 



J. 



(e) 



Thus the contribution of Eq. ( 4.30| ) reduces to 

\ M w) 3 2 



(4.32) 



(4.33) 



Despite appearances, the remaining terms in Eq. (4.26), which contain q in the denominator, 
yield a well-defined result in the (Q = 0, Q — > 0) limit. We relegate the details of the calculation 
to Appendix IB] and give the final result here, 



28^ 2 r (e) 

3 



(e) 



(4.34) 



including the contribution of the regular terms given in Eq. ( |4.33| ). 
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4.3.2 Computation of 

For G( B \ we can directly put q = and use Eq. ( |4.20|) in Eq. ( |3.17D . We then calculate 



1 



Tr^(0,0) 



d 4 p 



3 S(p 2 )rj e (p)l k 2 [k -pr/xp- 4k x p p - 4k p p x + 4p x p p ] 



2M^J (2vr) 
+ 4k-p [2k\p p + 2k p p x - p x pp - k x k p ] j . 

Using Eq. ( |4.19| ) now and carrying out the integral over p°, from Eq. ( |4.22| ) we obtain 



(4.35) 



G (B) (0,0) 



d 3 p 



8k ■ p 2k ■ v p ■ v — k ■ p (f e + f, 



J (2vr)32£ 
- 4k ■ p\(p ■ vf + (k ■ v) 2 } (f e -h) 



(4.36) 



This can be written as 



G {B) (0, 0) = ^- 1 16(k • v)k»v v Ip V - 8kWIp U - 4k^v x Ip uX - 4(fc • vfk^Ip \ , (4.37) 



and using Eq. ( 4.32| ) 



^(0,0) = j£- {-4^ + ^ 2 J^ - 4^} 



(4.38) 



4.3.3 Computation of G^ and G^ 

Proceeding in similar fashion, we get 



1 



2 Tr#Gg?(0,0) = JLf- J ^5(p 2 ) Ve (p)(p-k) 2 [ VX p(k-p)+p x kp + k xP p}. 



d 4 p 



(4.39) 



Using Eq. ( |4.19| ) and Eq. (4.22), we can now write 



G (c) (0,0) 



d?p 



M*J (2^f2E 



4k ■ p p ■ vk ■ v — k ■ p (/ e + f, 



4(k ■ v)k»v v Ip V - 4A^%„ , 



(4.40) 



and similarly it can be verified that G( D )(0,0) is given by the same expression. Using Eq. Q4.32P 
then 



G (c) (0,0) =G (D) (0,0) 



9 



Mfa} 3 



2 7 (e) 



(4.41) 
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5 Neutrino dispersion relation 



Collecting the results for the G (x) (0, Q -> 0) given in Eqs. Q4.34| ), (|4.38| ) and Q4.41] ) and substituting 
them in Eq. ( [4.21 ), we finally obtain 



4w 3 + 5H W 2 j(e) 

i 3 ^ 



icxt 



(5.1) 



for u = ±K, which is the result that must be substituted in Eq. ( [4.14| ) for the dispersion relation 
of the electron neutrino. For Eq. ( [4.14 ) holds but the term is absent. 

The integrals defined in Eq. ( 4.31 ) cannot be evaluated in a general way for arbitrary distribu- 
tions and rrif ^ 0, but for m e = they can be performed exactly to yield the following result in 
terms of the chemical potential and temperature of the electron gas: 



J 



(e) 



Hi p e T 



J. 



(e) 



12vr 2 

,4 



16vr 2 



+ 



12 ' 

M 2 T 2 77r 2 T 4 



+ 



240 



(5.2) 



(f) (f) 

For our present purposes, it is sufficient to note that, in any case, J\ and J2 are simply related 
to the number density rif and energy density pj of the particles in the background medium by the 
formulas 



-n 



if) 



j. 



!(„,-„,, 

1 



Thus, using Eq. ( |5.1| ), the dispersion relation becomes 



(5.3) 



(2) 



K + 2i^ ext (l + 7) + (1 + <f xt )& mat + b% 



K + 2K(j) ext (l + 7) + (1 + <f> ext )b 



l(2) 
} G ' 



(5.4) 



where 



7 



7 



4M^ 



26 

-2K(n e - Tie) + — (p e + Pe 

26 

2K(n e - Tie) + — (pe + Pe) 



for i/ P 



for v e , 



(5.5) 



with 7 (7) being zero for i/„ 
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The other terms in Eq. fl5,4| ) are the following. In Ref. |4|, we calculated the contributions to 

f2) 

b G for various kinds of background particles — non-relativistic or ultra-relativistic, non-degenerate 
or degenerate, arriving at the final formula given in Eq. (4.30) there. For the electrons, we use here 
the results for the ultra-relativistic case since it is the one that corresponds to the limit m e = 
that we have adopted in the present paper. Then, 

bg = -f^ <f xt x \ (5.6) 

for / = e,n,p, where Xf is the vector neutral-current coupling of the fermion, and 

J e = -5(n e - Tie) 

-f3mNriN classical nucleon gas 



"n^ degenerate nucleon gas , 



(5.7) 



FN 



for N = n,p, where vfn stands for the Fermi velocity of the nucleon gas. 

The term 6 ma t can be split into two parts. One is the usual Wolfenstein term, which is given by 



(n e - Tie) +Y,fXf{n f - nf) for v e , 
J2f Xf(n e - He) for v T 



~(2) (2) 

for the neutrinos, and fr^at = — ^mat f° r the antineutrinos. The second part consists of the correc- 
tions of 0(1/Myy). Using the results of Ref. ||] [in particular Eq. (3.27)], but putting m e — ► 0, and 
neglecting the terms that arise if there are neutrinos and other charged leptons in the background, 
we can write this part in the form 

/ _2 \ [ -\K(p e + Pe) foiV e , 



u ' I for V^Vr 



with foj^t = ^mat f° r the antineutrinos. 



6 Discussions and Conclusions 



As the formulas given in Eqs. ( |5.4D and ( |5.5| ) show, the effects of the induced gravitational terms 
on the neutrino dispersion relation are more noticeable at high energies. An important point to 
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remember is that the calculations that we have presented rely on the approximation made in Eq. 
@ for the VF-boson propagator. This implies that the results are valid as long as (p — k) 2 <C M^, 
which translates to K(E e ) <C Myy, where (E e ) denotes a typical average energy of the electrons in 
the background. For backgrounds with a temperature of the order of the electron mass or somewhat 
above, this implies that K <C W 7 GeV. Thus, if we consider values of K ~ 10 5 — 10 6 GeV , for which 
our results hold, we find that the induced gravitational effects can be substantial if the gravitational 
potential is not too small. In fact, from Eq. (|5.5| ), it follows in that case that 

for z^ e , with 7 = —7 for v e , and zero for the other neutrino and antineutrino types. Moreover, it 



is possible that, for a certain range of values of the neutrino momentum, the term in Eq. (6J) is 
the dominant one in the neutrino dispersion relation. These considerations may be applicable, for 
example, in the environments of Active Galactic Nuclei. It has been observed that in such systems, 
high energy neutrinos can have resonant spin-flavor transitions due to the combined effects of the 
gravitational interactions and the presence of large magnetic fields ||. Since in such environments 
0ext ^ 1/10 _ 1/100 in the region in which the neutrinos are produced, the momentum-dependent 
gravitational terms modify the neutrino dispersion relation in a substantial way, and therefore they 
can have observable consequences in the determination of the high energy neutrino fluxes from 
such systems. The calculations that we have presented, and our result for the neutrino dispersion 
relation, are useful in this and possibly other contexts, and provide a firm basis to study their 
possible effects in detail. 



Appendices 

A Proof of Eq. (^T5|) 

To prove the statement that (2v x v p — r] Xp )Y\ p (0, Q — > 0) is well defined, we consider the integral 
formula for Y\ p (£l, Q) given in Eq. ( [4.161) , Carrying out the integration over p°, in the rest frame 
of medium, and using the fact that, considered as a function of p and q, 



we have 



where now 



Yxp(n, Q) 



jx P (-p,q) = -j\ P (p,-q) , 
d 3 P 



(fe ~ fe) 



(2tt) 3 2£ 



3\ P 



q 2 — 2p ■ q 



+ (g^ -q) 



£ = \V\. 



(A.l) 



(A.2) 



(A.3) 
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We consider separately the functions 



Y 7 j 
Y 

1 vv 



v x vPY Xp , 



(A.4) 



and show explicitly that each one has a well defined value in the limit (O = 0, Q — ► 0), which proves 
the statement. 

From Eq. ( [4.17 ) we obtain 



il Xp 3\ P = -(.k-p)^(2p 



1 



(A.5) 



which implies that 



Y, 



d 3 P 
(2tt) 3 2£ 



(fe ~ fe)k-p 



-(n e - n s )k ■ v 



(A.6) 



For Y vv , setting D = 0we obtain 

i r d 3 P 



Y vv (0, Q) 



,., {2nf Ue-h){28^-28V-K) 



1 



1 



+ 



1 



(X)3 



(f e -f e )(SQ-K-<JP-Q) 



2V-Q-Q 2 2V-Q+Q 2 
1 1 



2V ■ Q - Q 2 2V-Q + Q 2 



.(A.7) 



The integrals with the factors 2£V ■ K and £Q ■ K are zero by symmetric integration since the 
corresponding integrands are odd under V — > —V. The remaining two integrals are evaluated by 
making the change of variables V —> V ± |<2 when the denominator is 2V ■ Q =p Q 2 , respectively, 
and then taking the limit Q — > 0. In this way, 



d 3 V 
(2vr) 3 



K/e - m 



(A.8) 



B The limit G^ A \U = 0, Q -> 0) 



The expression for ) appears in Eq. ( |4.26| ). In the text, we have evaluated the terms without 
any denominator. Let us denote the remaining terms by G^ A \ For these terms, we first put O = 0. 
The resulting expression for ^#(p, q) is written as 



^o(p, Q) 



-A£ s lo + A£ 2 V ■ K + 2£uV • Q - 2£ 2 K ■ Q . 



(B.l) 
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which follows directly from the definition in Eq. ( |4.27| ) . We can then write these potentially singular 
terms in the form 



& A \0, Q) 



9 



d 3 V 



(k 2 - 2k ■ p)^o(p, Q) f , £) _{k 2 + 2k 



(-P,Q) 



2V-Q-Q 2 



2V-Q+Q 2 



(k' 2 -2k' -p^o^-Q) 
2V-Q + Q 2 

G{ A) (0, Q) + G ( 2 A) (0, Q) 



fe(£) + 



(k' 2 + 2k' -p)^ (-p,-Q) 
2V-Q-Q 2 



(B.2) 



where G\ comes from the first pair of terms and G2 from the second pair. 

The integrals in this equation can be conveniently expressed in terms of the following ones: 



•F'i n,r 



as well as 



,{n) 



(2tt) 3 2£ 
d 3 V 



£ n (V ■ Q) r 



fe(£) 



(2vr) 3 2£ 



f n V- V ■ ■ ■ V 



2V-Q-Q 2 

fe(£) 



("I) 



n+m 



H£) 



2V-Q+Q 2 



2V-Q-Q 2 



\n+r 



f~e{£) 



2V-Q+Q 2 



(B.3) 



(B.4) 



where i\ etc are spatial indices. Using these notations, we can write the first pair of terms for G^ A > 



as 



g[ A) (0,Q) = 16cj 2 ^4,o + 8wi? • Q^3,o - 8w 2 ^ 2i i 



(B.5) 



where we have now used Eq. ( 4.1 9| ) . To simplify further, we first notice that by elementary argu- 
ments the integrals in Eq. ( |B.4[) can be expressed in terms of the integrals J- n ,m : 



(n) 


Qi 


i 


Q 2 


(n) 


1 " 




2 



F, 



n+2,0 



Q 2 



n+2,0 



3F ra , 2 \ QjQj 
Q 2 ) Q 2 



(B.6) 



The next step is to realize that the scalar integrals T n ^ m for a particular value of m can be expressed 
in terms of J- n - m for smaller values of m and the integrals Jn^ defined in Eq. ( [4.31 ). For example, 



4 e) + Q 2 ^,o 



(B.7) 
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Finally, in J- n p, the angular integration can be performed exactly to give 

1 



i ( l~{Q/28) \ 

lfiWo Q \l + {Q/2£)) 



MS) + (-i) n h{S) 



For small Q, the logarithm can be expanded in power series and we obtain 

r - 1 7- (e) Q2 f (e) +n(o*\ 

The integrals for larger values of m can be obtained by using Eq. flB.7| ). 

When we put the results of these integrals, Eq. (B.5) reduces to the form 



(B.8) 



(B.9) 



(B.10) 



In a similar fashion, the other pair of terms in can be evaluated. The result is 

Gi A \0,Q) = -12u^ - 4 £) + lM f 2 ' ~ 4 £) ~ ^4 £) + 0(Q 2 ) ■ (B.ll) 

Summing these two contributions and reinstating the overall factor, we can write 

-.2 



G (A) (0,0) 



9 



12u 2 4 e) + 4u;J 3 (e) 



(B.12) 



Adding this with the contribution of the regular terms obtained in Eq. ( |4.33| ) , we get the total 
result shown in Eq. (4.34). 

An alternative evaluation of Eq. (B.12) is as follows. In the Q-dependent integrals we make 
the change of variables V — > V ± in those terms containing in the denominator the factors 
V • Q T Q 2 ) respectively. By expanding the integrands in powers of Q and taking the limit Q — > 
we obtain in this fashion 



G (A) (0,Q^O) 



5" 



d 3 P 



M^J (2vr)3 2 f 



2Jfc 



P-i?| 2 ^(£/e(£)) 



d 

d£ 



(S - u) 2 £f e 



lv + K]2 ^s {£fs{£)) ~ lie l {£ + u)2£h 



d 3 V 



M^J (2vr)32£ 



4k ■ pi 



at 



k-p^-+S(S + u)f B 



(B.13) 
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using Eq. (4.19) in the last step. Carrying out an integration by parts and combining some terms, 
we get 



9 



d 3 V 

which can be written in the form 



I2(k ■ p) 2 (f e + h) + 4k ■ p{£ 2 - £u)(f e - h) 



g< a >(o,q-o) 



9 



Ylk»k v I^ v + 4k»v v v x I^ x - Auk^I^ 



(B.14) 



(B.15) 



in terms of the integrals defined in Eq. ( f4.29| ). Substituting into Eq. ( |B.15| ) the formulas given in 
Eq. ( p~3p , the result given in Eq. ( |B.12| ) is reproduced. 
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